A New 5 Flavour NLO Analysis and Parametrizations of 
Parton Distributions of the Real Photon 



F. Cornet", P.Jankowski 6 and M.Krawczyk 6 

a Departamento de Fisica Teorica y del Cosmos, Universidad de Granada, 
Campus de Fuente Nueva, E-18071, Granada, Spain 

b Institute of Theoretical Physics, Warsaw University, 
g ■ ul. Hoza 69, 00-681 Warsaw, Poland 

O 

Or 
< 

Abstract 

> . 

. New, radiatively generated, NLO quark (u, d, s, c,b) and gluon densities in 

a real, unpolarized photon are presented. We perform three global fits, based 
on the NLO DGLAP evolution equations for Q 2 > 1 GeV 2 , to all the available 
structure function F%(x, Q 2 ) data. As in our previous LO analysis we utilize two 
theoretical approaches. Two models, denoted as FFNScxftT & 2 NLO, adopt 
the so-called Fixed Flavour-Number Scheme for calculation of the heavy-quark 
contributions to i^(rE,(5 2 ), the CJK NLO model applies the ACOT(x) scheme. 
We examine the results of our fits by a comparison with the LEP data for the 
Q 2 dependence of the F^ , averaged over various x-regions, and the F^ c . Grid 
• parametrizations of the parton densities for all fits are provided. 
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1 Introduction 



In this work we extend the LO QCD analysis of the structure function F^(x, Q 2 ) of the 
unpolarized real photon [1-3], to the Next-to-Leading order (NLO). As before we are 
especially interested in the description of the charm- and beauty-quark contributions to 
F^(x, Q 2 ). We adopt two approaches, the FFNScxr: and CJK ones, as discussed in [1] 
and (with some modifications) in [2]. In the first model, referred to as FFNScj^I 
NLO, we adopt a standard approach in which a heavy quark, h, contributes to the 
photon structure function only through the 'direct' (Bethe-Heitler) 7* + 7 — » h + h 
process and the process with 'resolved photon' 7* + G 7 — > h + h. In the FFNSc je- 2 
NLO model we include in addition 7* + G 7 — » h + h + G and 7* + g 7 (<f r ) — > q(q) + h + h 
processes 1 . In the third model, CJK NLO, we use the ACOT(x) scheme [4] in which, 
apart from the direct and resolved photon 0(a s ) contributions 2 , we deal with the 
heavy-quark densities qh{x,Q 2 ). The subtraction of the double counted terms as well 
as introduction of a new kinematic variable Xh improve description of F^(x, Q 2 ) in 
both regions: below and above the heavy-quark thresholds. 

In [1-3] we followed the idea of the radiatively generated parton distributions in- 
troduced by the GRV group first to describe nucleon [5] and pion [6] and later used for 
the real [7], [8] and virtual photon [8], [9]. In our LO analyses [1-3] we used Qq = 0.25 
GeV 2 and adopted the input distributions of the valence-type form introduced in [7] 
neglecting sea-quark densities. In the case of the NLO analysis we found that the 
starting scale of the evolution appears to tend to higher values, Ql ~ 0.7 GeV 2 (simi- 
larly as in [10]). We tested the assumption on vanishing sea-quark densities at Ql by 
performing additional fits. 

All our global fits are performed at Q 2 > 1 GeV 2 utilizing the set of all available 
F 2 7 (x,Q 2 ) data. 

We test our results for F^{x, Q 2 ) by comparing them with the LEP data which were 
not used in the fits: for the Q 2 dependence of the F^ , averaged over various x-regions, 
and for the data. 

Our paper is divided into six parts. In section 2 we shortly recall basic schemes 
applied to calculation of the heavy-quark production in the Deep Inelastic Scattering 
on the real-photon target. Section 3 is devoted to the description of the FFNSc jk NLO 
models of the F 2 7 (x, Q 2 ). Next, the CJK NLO model based on the ACOT(x) scheme [4] 
is presented in detail. In section 4 we describe the solutions of the DGLAP evolution 
equations in all models and assumptions for the input parton densities. Results of the 
global fits are discussed and compared with the F^x.Q 2 ) , F^iQ 2 ) and F^ data in 
the fifth section of the paper. We end our paper by presenting in the Appendix the 
technical details of the calculation. The parton distributions resulting from our fits 
have been parametrized on the grid. 

1 These two FFNSc jkI & 2 models differ from the corresponding models introduced in our LO 
analysis. 

inclusion of the 0(a 2 s ) terms in the CJK NLO approach is beyond the scope of this work. 
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Figure 1: Deep Inelastic Scattering on quasi real photon (P 2 ~ 0), e + 7 ^ e + hadrons. 

2 Production of heavy quarks in DIS e7 

The Deep Inelastic Scattering on a photon (DIS e7 ) allows to measure the structure 
function K7, and others, via the process 

e + 7 — > e + hadrons, (1) 

presented in the figure H The differential cross-section for this process is given by: 



dxdy xyQ x 



[(1 + (1 - y) 2 )F2(x, Q 2 ) - y 2 F2(x, Q% (2) 



where x and y are the standard DIS e7 variables, defined through the virtuality Q of 
the probing photon and momenta of interacting particles shown in figure d 

x - v - ^ (3) 

2p ■ q p ■ k 

and F2{x,Q 2 ) is the so-called longitudinal structure function corresponding to the 
exchange of the longitudinally polarized 7*. For the presently kinematically accessible 
values of y the F2(x, Q 2 ) contribution is negligible. 

Each of partons (quarks and gluon) of the real photon contributes to the F^(x, Q 2 ) 
proportionally to the corresponding partonic cross- sect ion. The calculation of the light 
u-, d-, s-quark and gluon cross-sections is straightforward. 

Various schemes used in the calculation of the heavy-quark production have been 
described in [1]. They are also explained in the vast literature, [4], [11] and [12]. Let 
us shortly recall the main facts on the FFNS and ACOT(x) approaches applied in this 
analysis. 

There exist two standard schemes in a calculation of a heavy-quark, h, production 
at the hard scale fi. In both of them the light quarks u, d and s are treated as massless 
because for them /i ^> Aqcd > m q , where m q is a light-quark mass. In the case of the 
DIS process the hard scale \i is identified with the virtuality of the probing photon, 
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/i 2 = Q 2 . In the massive Fixed Flavour-Number Scheme (FFNS) heavy c- and b- 
quarks are treated differently: the massive charm and beauty quarks produced in hard 
subprocesses 7* + 7 — > h + h can only appear in the final state. In the second, massless 
scheme, called the Zero-mass Variable Flavour-Number Scheme (ZVFNS) heavy-quark 
distributions appear alike the light-parton densities. When Q 2 is larger than a threshold 
associated with a heavy quark (usually taken as Q 2 = ml), this quark is considered as 
an extra massless parton in addition to the three light quarks. This way, the number 
of different types of quarks (flavours) increases with increasing Q 2 . One introduces a 
notion of "active quarks" , for which the condition Q 2 > m 2 is fulfilled and which can 
be treated as massless partons of the probed real photon. When active, their densities 
qi(x,Q 2 ) differ from zero, otherwise qi(x,Q 2 ) = 0. 

The above standard schemes are considered to be reliable in different Q 2 regions. 
The FFNS looses its descriptive power when Q 2 ^> m 2 , on the other hand the ZVFNS 
does not seem appropriate if Q 2 ~ m 2 . To obtain a prescription working for all hard 
scales one needs to combine features of both schemes. There exist the whole set of 
the approches, generally denoted as General-mass Variable Flavour-Number Schemes 
(GVFNS), which in the Q 2 — > m\ and Q 2 — > 00 limits reproduce the behaviour of 
the FFNS and ZVFNS schemes respectively. Their recent realization is the so-called 
ACOT(x) scheme introduced for the proton structure function in [4] and used by us 
in our LO analysis [1-3]. 

In the ACOT(x) scheme, all contributions to the considered process, which would be 
included separately in the ZVFNS and FFNS, are taken into account. Such procedure 
requires proper subtraction of double counted contributions, i.e. the large logarithms 
which appear in the massive heavy-quark Wilson coefficients (the FFNS scheme) but 
are also resummed in the qh{x, n 2 ) densities. That way one obtains a scheme in principle 
appropriate in the whole Q 2 range. Still, another problem emerges. In any process the 
kinematical threshold for the heavy-quark production is given by the total energy 
of that process, W. Obviously W must be greater than two masses of the heavy 
quark in hand, W > 2m h , so in the DIS e7 case where the hard scale is Q 2 , we have 
W 2 = (1 - x)Q 2 /x > 4m);. As /j 2 = Q 2 the ZVFNS condition \i > m h (or any 
similar) for treating a heavy quark h as a massless parton is not correct. It may 
happen (for small enough x) that qh{x,Q 2 ) = in the kinematically allowed region 
W > 2mj L . On the other hand, we can also obtain nonzero heavy-quark densities in 
the kinematically forbidden (x, Q 2 ) region. The ACOT(x) scheme solves this problems 
through introduction of a new variable Xh = ^(1 + ^ m h/Q 2 ) which replaces in the 
ZVFNS the Bjorken x as an argument of the heavy-parton density. 

In the next sections we will present the realizations of the FFNS and ACOT(x) 
schemes for the calculation of the photon structure F£ in the Next-to-Leading order of 
QCD, namely our FFNS CJif NLO and CJK NLO models. 

3 Description of the F 2 7 in the FFNS CJA NLO 

This section is devoted to the FFNScxrT & 2 NLO models, which are the realizations of 
the FFNS approach for the photon structure F^ in the Next-to-Leading Order of QCD. 
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The only difference between this two models is following: the FFNScj^NLO 2 approach 
includes additional higher order contributions to neglected in the FFNScjif NLO 1 
model. 

First, we describe the light-quark contributions to F^x^Q 2 ) which are common 
for the both analysed models. Next, we discuss the heavy-quark contributions for the 
FFNScji^NLO approach. Further, we recall the DIS 7 factorization scheme used in the 
model and resume by giving the final formulea for the F^ix, Q 2 ) in the FFNS^xft- NLO 
models. 



3.1 Light-quark contributions 

In the Next-to-Leading Order logarithmic approximation of QCD or, in short, in NLO 
QCD the light-quark contributions to the photon structure function F^x^Q 2 ) can be 
written in terms of light-quark (u, d and s) densities q](x, Q 2 ) as follows: 



iitfOr, Q%i 9h t = £ e 2 {(q] + %)(x, Q 2 ) + e^C^x) 

i=l 



(4) 



a s {Q 2 
2tt 



V 



(Ql + Ql)(y,Q 2 )C^(-) + G^ yi Q 2 )Ci%(-] 



Note, that q](x,Q 2 ) = qJ(x,Q 2 ). The C^Ax) functions are the 0(a\) order terms of 
the hadronic (Wilson) coefficient functions (j = 7, q, G) 



C 2d (x, Q 2 ) = Ci%) + ^^C^ix) + -... 



(5) 



Each of them is related to the hard process of a parton originating from the real photon 
(for instance quark for C 2 , ? (^)) interacting with the virtual photon. The coefficients 
appearing in Eq. (j3J) are listed in table ^ The formulae for the C^j (x) functions up 
to the 0(a 2 ) order can be found in [13], see also the Appendix. 



order 


hard (parton) process 


coefficient function 


O(a s ) 




C^(x) = 6(l-x) 


7* + 7 — ► q + q 


C^(x) 


O(al) 


l* + q^m^q(q) + G 






C 2G (x) 



Table 1: Hard parton processes and the corresponding Wilson coefficient functions for light 
quarks. 



The evolution of the parton densities in InQ 2 is governed by the inhomogeneous 
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DGLAP equations [14]. In the NLO we have: 



dq](x ' Q) = ^k,{*,Q 2 ) 



+ 



dlnQ 2 
a s (Q 2 ) [Uy 



+ 



2tt 

dG^(x,Q 2 
dlnQ 2 

a s (Q 2 ) [Uy 



N, 

E 

fe=i 



ql( y , Q 2 )P qg (^ Q 2 ) + ? k (y, Q 2 )P q «(-, Q 2 " 



Ot 

—k G (x,Q 2 ) 

ZTX 



2tt 



Nt 



fe=i 



?2(i/, Q 2 )^ 9 (~, g 2 ) + fib, Q 2 )p G ,{^ q 2 ) 



(6) 



. x 



G^(y,Q')P qG (-,Q 2 

y 



(7) 



G\y,Q 2 )P GG {-,Q 2 ) 
y 



where Nj is a number of active quarks. Here Nf = 3 as all heavy-quark densities are 
equal zero in the FFNS models. The functions Pij(x,Q 2 ) and ki(x,Q 2 ) are the NLO 
splitting functions 



h{x,Q 2 ) = kf\x) + 



a s (Q 2 
2tt 



kP(x) 



(8) 



P l3 (x,Q 2 ) = P$>{x) 



^pp. ( M 

Ztx j 



The k q (x, Q 2 ) and k G (x, Q 2 ) are reffered to as the photon-quark and photon- gluon 
splitting functions, respectively. Up to the NLO order P Gq = P Gq . Formulae for 0(a s ) 
and 0(a 2 ) splitting functions can be found for instance in [15]. 



3.2 Heavy-quark contributions 

In the FFNS models heavy quarks appear only in the final state of the process and the 
density functions qh{x,Q 2 ) are equal to zero. The charm- and beauty-quark contribu- 
tions to the structure function are described by the heavy-quark coefficient functions 
C*2 j related to the hard processes as listed in table 121 One obtains 



-F2{x,Q 2 ) 



\heavy,FFNS 



E 



' It 



h(=c,b) 

„i , r s 



a s (Q 2 ) [Uy 



2tt 



/ d J>G^y, Q*)CZJ\ X -, %) + «j[f C*f (*, % 
J Xh y • V rn{ 2w m' h 



E 

h(=c,b) 

y rn h . 



2?r / J X h y ~t v m h y m h 



)) 



| aa s {Q 2 ) c 4 C ft,(i) (x ^ 



(2tt)* 



mi 



where the parameter Xh = ^(1 + ^ m h/Q 2 ) takes care of the proper region of the 
integration. 

We consider two FFNS models, FFNScj^I and FFNS CJif 2 NLO, where only in 
the latter we include the 0(a 2 ) and 0(aa s ) contributions to F%. The reason for 



(9) 
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order 


parton subprocess 


coefficient function 




Y +1 ^ h + h 


/Vl,(0) 
U 2,7 


0(al) 


7* + -> h + h 


i, / 1 \ 
°2,G 


Y + 1 ^ h + h + G 






Y + G J -*h + h + G 


2 G 


7* + <f (<f) -> q(q) + h + h 


/AW ^(2) 
°2,</ ' U 2,„ 



Table 2: Hard parton processes and the corresponding Wilson coefficient functions for heavy 
quarks. 



considering two models is twofold. First, to check the size of the higher order heavy- 
quark contributions. Secondly, to perform a test whether the 0(af) and 0(aa s ) terms 
should be included in the calculations as they are of higher order than the other NLO 
contributions. 

Let us notice that among the higer order terms there appear two Wilson coefficients 
originating from the process 7* + q 1 {(f) — * q{q) + h + h, the C^'g and (see Fig. |2J), 
the interference term of two diagrams (proportional to e^q) gives according to [13]. 
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h 





Figure 2: Diagrams relevant for 7* + q' 1 {q l ) — > q(q) + h + h process corresponding to a) the 



C^'i coefficient and b) the C^l coefficient. 



(2) 



"2,g 



3.3 The MS and DIS 7 factorization schemes 

All the splitting and coefficient functions introduced in former sections have been de- 
rived in the modified minimal subtraction factorization scheme (MS). (They all can 
be found in the literature, see Appendix A). Let us notice that the C^hx) term in Eq. 
IH describing the 'direct' contribution of the photon to F^x, Q 2 ) causes a problem in 
the MS scheme [16]. Firstly, contrary to all other contributions, this term does not 
vanish at the input scale Qq. One obtains very different boundary conditions for the 
light-quark contribution in LO and NLO: 

-/^Ql)\ Ught ={ E i ie l c m {x) ;nlo. (10) 
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This leads to very distinct results over the whole x-range (see for instance [17]). The 
heavy-quark contributions, C^^(x, £2-), for k = 0, 1, of Eq. (jHJ) do not cause a similar 

problem as the 0(a®) term appears already in the LO calculation and the C^'^fa, ^-) 
coefficient is only its 0{a\) correction. 

Secondly, the C2 (%) term calculated in the MS scheme is negative and divergent 

in the large-re region, where C^(x) ~ 3[ln(l — x) — 1], leading to unphysical negative 
values of K7 at x — > 1. It is not cured by addition of other NLO contributions to the 
structure function as they all vanish at x — 1 (if a typical choice of the shape of the 
hadronic input is made). 

In a new factorization scheme, denoted as DIS 7 , introduced in [16], one avoids the 
troubles connected with the Cf^x) term by absorbing it into the quark distributions 
q]{x, Q 2 ), namely 

?r^,Q 2 ) = <ir s M 2 )+e^c^ m (x), (ii) 



C^(x) = ^' MS (*) + A^(*) = 0, 



where obviously 



ACi%) = C^ MS (x). (12) 

The gluon density and other coefficient functions are unaltered by the above redefinition 
but it leads to a modification of the photon-quark and photon-gluon splitting functions 
appearing in the DGLAP evolution of the parton distributions: 

k^(x,Q 2 ) = ^MV^^ (13) 

# dis ^,q 2 ) = # MS (*,g 2 )-E e < / d M°l(-)c^(y)- 

In our analysis we apply the DIS 7 factorization scheme, for simplicity the DIS 7 
superscript will be omitted below, if possible. 

3.4 Formulae for F 2 7 

F^(x,Q 2 ) in the FFNS CJK NLO 1 model reads; 
1 



F2(x,Q 2 )\ FFNS l= (14) 
x 



dy 

y 



s 



In the case of the F%{x, Q 2 )\ffns^- model we include the 0(a 2 ) and 0(aa s ) terms, 
obtaining 

_ 1 



Fl(x,Q 2 )\ FFNS 2 = -F2(x,Q 2 )\ FFNS l 
x x 



(15) 



+ £ 

h(=c,b) 



2tt 



2\\ 2 /•! 



- e ^(,,g 2 )^(^. 



2 



£(?7 + Q 2 )(e^£ } (-, ^) + efcjf (*, ^)) 

,y y m h y m h 



i=l 



aa s (Q 2 ) 4 /!,(!), Q' 2 



{2k) 



ni, 



4 Description of the in the CJK NLO model 

This section is devoted to the description of the CJK NLO model based on the 
ACOT(x) approach. The approach combines the FFNS scheme introduced in the for- 
mer section with the ZVFNS scheme in which heavy quarks are treated in the similar 
way as the light ones. 



4.1 Light-quark contributions 

The light-quark contributions to F% are the same in FFNS and ZVFNS schemes and 
were already given in Eq. (j3J) for the MS factorization scheme and in Eq. ()14|) for 
the DIS 7 scheme. The difference between the numerical values of those contributions 
originates from the difference between the number of the active quarks, Nf, in the 
DGLAP equations in both approaches. We have Nf = 3 for the FFNS^xft- NLO 
models, while Nf = 5 in the CJK NLO model. 



4.2 Heavy-quark contributions 

Apart from the heavy-quark terms which properly describe the /i-quark contributions to 
F^(x, Q 2 ) , in the corresponding Q 2 w ml regions, are given by the FFNS formula (J9*|). 
Let us only notice that in the CJK NLO model case we neglect the 0(a 2 ) and 0{aa s ) 
contributions to F%{x, Q 2 )\heav V ,FFNS- They lead to a very complicated expression for 
F^x, Q 2 ) and their analysis is beyond the scope of this work. As already discussed, the 
FFNS scheme looses its descriptive power if Q 2 ^> m 2 h . Therefore, the addition of the 
ZVFNS terms, the massless charm- and beauty-quark distributions, qh(x,Q 2 ), in the 
photon, appropriate in that region, is needed. They are combined with the light-quark 
densities through the DGLAP equations with Nf = 5 active quarks. In the ZVFNS 
approach heavy-quark contributions to F^x.Q 2 ) read 



-F](x,Q 2 



+ 



2vr 



'\heavy, ZVFNS 

l dy 

y 



z 

E 4{(qI + €)(^Q 2 ) + 4^%) (i6) 



h(=c,b) 



(<ll + €)(y,Q 2 )C^(-) + G^(y,Q 2 )C^) 
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4.3 Subtraction terms 



A simple summing of all the discussed terms, namely 

F^(x,Q 2 ) = F^ix, Q 2 )\light + F^ix, Q 2 )\heavy,FFNS + F^(x, Q 2 )\heavy, ZVFNS, (17) 

would lead to double counting of some contributions. Let us notice first that the ZVFNS 
subprocess with the direct coupling of the real photon to the (massless) heavy quark, 
term C^(x) in Eq. (fTEj). was also included in the FFNS formula as C^j (x, 

Those two terms originate from the same 7* + 7 — > h + h process, however they are 
not equal as the c!Q(x) term is obtained in a massless approximation. The same 

is true for the Cf^x) (Eq. (fTS|)) and C^'q (x, ^-) (Eq. 0) terms arising from the 

7* + G 7 — > process. Obviously, only one set of the above terms should be included 
in the final CJK NLO formula for the F^x^Q 2 ). We choose to apply the massive 
C2^°\x, -^t) and C2q\x, formulea, as unlike the massless approximations those 
terms properly vanish at the heavy-quark threshold. 

The C^g coefficient in Eq. is related with the process 7* + h? — > h + G, with a 
heavy quark in the initial state and therefore appears only in the ZVFNS approach. It is 
the only heavy-quark Wilson coefficient which we apply in the massless approximation. 

The massless C^hx) and C^q(x) terms are not the only ones that must be sub- 
tracted. Still, we double count the corresponding collinear configurations which are 
a part of the DGLAP equations and manifest themselves via the splitting functions 
Pij and ki. Such configurations give rise to terms proportional to InQ 2 and through 
DGLAP equations are part of the heavy-quark densities q^ix, Q 2 ). The same collinear 
configurations appear in massive heavy-quark Wilson coefficients C^j^. Therefore, we 
must subtract the overlapping contributions. In the case of our model, where we ne- 
glect the higher order heavy-quark Wilson coefficients, there are only two such terms, 
same as in the LO case described in [1] and [2] . They are related to the 7* + 7 — > h + h 
and 7* + G — > h + h processes. We calculate the necessary subtraction terms from the 
approximated integration of the corresponding parts of the DGLAP equations (J7J) and 
JSJ) over the region (Ql,Q 2 ) (see [2]). The integration produces the In §2 • e 2 h ^k q °\x) 

and In • ^ip £, 7 P${~ )G(y, Q 2 ) terms 3 . These terms muliplied by 2x and e 2 h , 
according to the formula F%(x, Q 2 ) = J2 q xef^q 1 + cp)(x, Q 2 ), must be subtracted from 
the sum of Eq. (JT7j). 

The described summation and subtraction procedure leads to the following formula 

3 Notice that we neglected the Q 2 dependence of the strong coupling a s and the gluon density 
G(y,Q 2 ). 
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Figure 3: Graphical representation of the contributions to the CJK NLO model. The first 
diagram represents the ZVFNS contribution. The third and fifth diagrams show the FFNS 
resolved-photon and direct-photon contributions, respectively. The second and forth dia- 
grams are the corresponding subtraction terms. 



for the F 2 7 (x,Q 2 ) : 

-F^(x, Q 2 ) = J2 {(?7 + Q 2 ) + e?fcS 



(0), 

2,7 



i=l 



+ 



cqg s 

2tt 



y 



+ 



E el\{ql + f h )(x, g 2 ) + e^Cjf (*, 

Z" 1 —G 7 (y, Q 2 )Pga(-)\ ■ 

Jx y 9 y J 



h(=c,b) 

a s (Q 2 
2vr 



- lng-2 e 2 ^)(x)-lng.2^^ 



gg 



gg 



Its graphical representation is given in Fig. El If one includes the 0(a 2 ) and 0(aa s ) 
heavy-quark contributions, the structure of the subtraction terms becomes more com- 
plicated, since one has to subtract terms proportional to k^\x), P^(x) and Pqq(x). 
Moreover these functions are divergent at x — 1 and therefore their integrals lead to 
numerical instabilities. Solution of that problem is beyond the scope of this work. 



4.4 The MS and DIS 7 factorization schemes 

The introduction of the heavy-quark contributions to the F^x, Q 2 ) forces us to apply 
a modified DIS 7 scheme in the CJK NLO model case. Let us recall that in this model 
apart from the massless, divergent at x — > 1 contribution of light-quarks, C^hx), the 
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massive, finite C 



2,7 1 



term appears as the analogous heavy-quark contribution. 



In order to remove the large-x divergence of the C. 



(0) 
2,7 



we proceed in the same way 



as described in the FFNS approach, absorbing that term into the light-quark distribu- 
tion functions q](x,Q 2 ). What about the massive C^^ix, £2-) contribution and the 
heavy-quark densities? The charm- and beauty-quark distributions in the CJK NLO 
model appear as a part of the massless ZVFNS scheme. The evolution of all partons is 
performed through the common set of five-flavour DGLAP evolution equations. There- 
fore, the subtraction of the C^^x) term from the Eq. (|18|) and modification of the 
k qi and kc splitting functions (expressed in Eq. ()13J)) affect the heavy-quark densities 
in the same way as they affect the light-quark distributions. Namely, the following 
redefinition occurs: 



q h 7 (x,Q 



^.(O^DIS-w ^ 



2 a ^(0),MS, 



U 2, 7 1 



Q 



m 



(19) 



As ACg(i) 



-a 



(0),MS 
2,7 ' 



X) 



see Eq. (JT2J), the last equality takes form 



^,(0),DIS 7 ^ = ^,(o),ms^ _ g (o),ms^ 



Q 2 



mi 



"> 2 . 
mi 



'2,7 



(20) 



In practice the subtraction of Eq. (|2Uj) can not be performed in such a simple way. 

Again, the C^' MS (x) term on the right side would lead to a divergence at large x and 
the advantage of the DIS 7 scheme would be lost. We have to proceed in a different 
way. 

We decide to resolve the above problem by loosing part of the information brought 
by the massive MS (x) contribution. That term calculated in the m\ « approx- 
imation leads to (see Appendix A) 



a 



,(0),MS/ Q 



2,7 



X, 



nil 



+ 



C^ MS (x) + et2k^ MS (x)\n 



mi 



Set 



-^(l-x)^+(x(l-3x) Aml 



Q 2 



Q- 



— X 



5m; 



(21) 



In 



1 + 
1-/3 



where we keep the heavy-quark mass and (5 



us notice now that the term ~ kq MS { 
term in Eq. (fTSj) only by the constant 
the MS scheme: 



Ami x 



X) 



In 



wherever it is possible. Let 



(l-x)Q 2 

of Eq. (|21|) differs from the first subtraction 



We omit that constant and obtain in 
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2." 



■r. 



-In 



mi 



Ql 



C { T\x) + Zei 
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which according to Eq. (j2*Uj) leads to the following relation in the DIS 7 scheme 



2n 2 ^ {X 'ml } lU Ql ^2^ {X) 



(23) 



IT 



_ , ,4m? / . ,4m? 



x 



imi 



Q 4 



In 



1 + 
1-/3 



In the following we use only the DIS 7 factorization scheme, and omit the superscript 
for simplicity. 



4.5 The Xh variables 

Let us recall that VFNS schemes widely utilize an approach in which ql(x, Q 2 ) = for 
all Q 2 values smaller than ml and q^(x,Q 2 ) becomes nonzero for Q 2 > m\. At each 
such step the number of active quarks grows by one. This approach leads to unphysical 
heavy-quark densities near their production thresholds. The kinematical threshold of 
the heavy-quark production, and therefore its non-zero density in the photon, is given 
in DIS e7 by the total energy of the ej* process, W. As discussed in section 2 W must 
be greater than the sum of two masses of the heavy quark in hand, W > 2m/,. W and 
Q 2 are connected, W 2 = (1 — x)Q 2 /x, which is the reason why the ZVFNS condition 
Q 2 > wn\ for treating the heavy quark ft, as a parton is too simple. 

In the ACOT(x) scheme used by the CJK NLO model the number of active quarks 
is equal Nf = 5 at the whole range of the DGLAP evolution. To ensure that a 
heavy-quark, ql(x,Q 2 ), distribution disappears when W — > 2m^ one introduces a new 
variable Xh = x(l + 4m 2 /Q 2 ) and calculates ql(XhiQ 2 ) m place of the q1(x,Q 2 ). At 
the threshold, when W = 2mh, we obtain Xh — 1 an d the corresponding heavy-quark 
density vanishes as desired. The form of the Xh variables could be chosen differently if 
only their above threshold behaviour was preserved. We chose to apply the same form 
as appears in the integration of the massive heavy-quark Wilson coefficients with the 
parton distributions, see also [4]. That way the same variable takes care of the proper 
vanishing of the /i-quark density and of the ~ ^■G 1 {y, Q^C^'q (| , ^r) contribution 
to F^x^Q 2 ) at the kinematic limit. For the consistency of the approach also in the 
subtraction ~ y-P g g '(^)Gr 7 (y, Q 2 ) term we exchange x with Xh and that way enforce 
its proper kinematic behaviour. 

As described in [1] and [2] the subtraction ~ k\ (x) term does not vanish as 
W — > 2rrih even when one replaces x with Xh- In the NLO case we avoid that problem 
because, as discussed in section 14 A\ that contribution cancels in the modified DIS 7 
factorization scheme. 

Unfortunately, a particular heavy-quark contribution appears in the NLO anal- 
ysis which causes new problems at the threshold. The contribution ~ J x -^(ql + 

Oh) (lli Q 2 )^2 q !(~) from the ZVFNS scheme could be rewritten in two ways, as ~ 
IL + €)kQ 2 )CS(§) or as ~ £ d f(ql + #)( 7fc , Q 2 )C^ y ), with lh defined 

analogically to Xh (lh = 2/(1 + ^ m h/Q 2 ))- The first approach leads to numerical insta- 
bilities for y values close to x originating from the ^jz^r In term in C^l (|). Second 
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form can not be adopted for the technical reasons in the Mellin space. Therefore, we use 
the unmodified formula of Eq. (|18|) . Finally, to avoid the non-zero heavy-quark contri- 
butions in the kinematically forbidden (x, Q 2 ) region we impose by hand the following 
constraint: ^-(qh + <jh)(y, Q 2 )C^q(^ ) = for Xh > 1- Obviously, the same reasoning 
would apply to higher order coefficients neglected in the CJK NLO model. 



4.6 Formula for 

Summing the light- and heavy-quark contributions, converting them into the DIS 7 
factorization scheme and finally introducing the Xh variables we obtain the final formula 
of the F 2 7 (x, Q 2 ) in the CJK NLO model: 



-Fl(x,Q 2 
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(24) 
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Q 2 
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1-/5 



To prevent the unphysical situation, in which the ZVFNS scheme contribution to 
F^txiQ 2 ) is negative, or in other words in which the total F^(x, Q 2 ) is smaller than 
the FFNS scheme contribution, as in [1-3], seperately for each heavy quark, we impose 
an additional positivity condition, here in form: 



(ql+qI)(x h ,Q 2 > 



'2tt 



iql + ql)(y,Q 2 )a 



2, 9 ' 



1 9 



^(,,Q 2 )P g ( 2(- 

Xh y y 

(25) 



> 0. 



5 Global fits - solving the DGLAP evolution 

In this section the technical details of the solution of the DGLAP evoultion equations 
are shortly described. 

5.1 Mellin moments 

The DGLAP evolution equations (J7J) and (JSJ) as well as the structure function formulae 
(|14J) and (}24|) are very much simplified if they are transformed into the Mellin-moments 
space. The n-th moment, where n is a complex number, for any function /, gener- 
ally depending on x and Q 2 is defined by f n (Q 2 ) = J x n ~ l f{x,Q 2 )dx. We aplly 
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this definition to the parton distributions, splitting functions and the light-quark Wil- 
son coefficients C^}{x). Their convolution integrals in the (x, Q 2 ) space transform 
into multiplication of the corresponding moments in the Mellin space. This simple 
correspondence does not hold in the case of the convolutions with the heavy-quark 
Wilson coefficient C^'q (x, 0?) and the P^q(x) splitting function because of the inte- 
gration limit Xh- Therefore those integrations must be performed in the (x, Q 2 ) space. 
All results calculated in the Mellin space are transformed through the inverse Mellin 
transformation to obtain their values in the (x, Q 2 ) space. 



5.2 Nonsinglet and singlet parton densities 

In the photon case we define the singlet (S 7 ) and non-singlet (/ns(jv/)) °L uar k distribu- 
tions as 



X 7 (x,Q 5 



1=1 

N f 



i=i 



with 



(26) 
(27) 

(28) 



i=l 



For the singlet density we sum over all quark flavours appearing in the evolution, the 
active quarks, Nf = 5. Although in CJK model Nf = 5, for the non-singlet case 
it is necessary to calculate all, <2ns(3)> ^ns (4) ano - ^NS(5)( a; ' *° obtain all six parton 
distributions, see [1]. 

The evolution of the non-singlet distributions is governed by the simplified DGLAP 
equation. In the Mellin moments space it reads 



dP' n 



NS(Vy 



,(Q 2 



a in _| a s(Q 2 ) p NS,n jr/,n 

27r K NS(iV / )+ 2yr r qq /nS(JV,) 



(Q 2 )- 



(29) 



dkiQ 2 

The singlet distribution evolution is coupled with the evolution of the gluon density 



d^> n {Q 2 ] 
dlnQ 2 

dGn> n {Q 2 ) 



2tt s 2tt 
2tt g 2tt 



dlnQ 2 

Thererefore, one has to solve the matrix equation 



P^ n (Q 2 ) + P£ G G^ n (Q 2 ) 
P^ n (Q 2 ) + P£ G G^ n (Q 2 ) 



(30) 
(31) 



dlnQ 2 



^kl 
2n - 



2tt 



(32) 
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with 



-fy.n 



h v 



pn pn 
EE EG 

pn pn 

GE ^GG 



The new point-like splitting functions are defined as 
k mNf) (x) = 2N f (<e 4 ) - (e 2 ) 2 ) 
kz(x) = 2N f (e 2 ) 



fc( o) , a stQ 2 ) h m 

q + 2vr 9 



fr (0) , ) ;/!) 



[X 



(33) 



(34) 



The formulae for the above spliting functions Pij(x,Q 2 ) can be found for instance 
in [15]. Their Mellin moments were calculated in [18]. 



5.3 Point- and hadron-like parts 

The solution to the DGLAP equations can be divided into the so called point-like (pi) 
part, related to a special solution of the full inhomogenous equations and hadron-like 
(had) part, being a general solution of the homogenous equations. Their sum gives the 
parton density in the photon, so we have 



/NS(iV / ),had( ( 5 ) + fm(N f ),pl(Q 

flUQ 2 ) + &Q 2 )- 



f~NS(N f )(Q 

fs \Q J — ^S,had 

The non-singlet hadron-like solution is given as 

,(0),n 



and the non-singlet point-like solution is 
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(l),n _ A , . (Q),w 
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(35) 
(36) 



(37) 



(3* 



(39) 



(40) 



Here /3q and /?i are the coefficients of the running strong coupling, a s , which we discuss 
in detail in section IB~5l 

The singlet point- and hadron-like solutions have the same form but are more 
complicated than the corresponding non-singlet solutions. All formulae can be found 
for instance in [17]. 
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5.4 Input parton densities. VMD 

The hadron-like parts of the singlet and non-singlet densities, denoted in this section as 
/had, need input distributions. For this purpose we utilize the Vector Meson Dominance 
(VMD) model [19], where 

/L(*>Q2) = £^/ V (*,QS)> ( 41 ) 

v J V 

with sum running over all light vector mesons (V) into which the photon can fluctuate. 
The parameters f v can be extracted from the experimental data on T(V —+ e + e~) 
width. In practice we take into account the p° meson while the contributions from the 
other mesons are accounted for via a parameter k 

J P 

which is left as a free parameter in the fits. 

We use the input densities of the p° meson at low Ql i n the form of valence- 
like distributions both for the (light) quark (v p ) and gluon (G p ) densities. All sea- 
quark distributions (denoted by ( p ) are neglected at the input scale. At this scale, the 
densities v p , G p and ( p are related, according to Eq. (|4^|) to the corresponding densities 
for a photon, see below. 

The v p density is given by 

v p (x,Ql) = \{u p+ +u p - +d p - +d p+ )(x,Ql), (43) 

where from the isospin symmetry 

u p+ (X, Ql) = u p - Or, Q%) = dT (x,Q 2 ) = d p+ (x, Qg) . (44) 

Note that all the densities in Eq. f)43j) are normalized to 1, eg. u p+ dx = 1. 

Two constraints should hold for the v p (x, Ql) density, the first one is related to a 
number of valence quarks in the p° meson, 

l 

2v p (x,Q 2 )dx = 2, (45) 
the second constraint represents the energy-momentum sum rule: 

/ x (2v p (x, Ql) + G p (x, Ql)) dx = 1. (46) 
Jo 

We parametrize the input densities as follows: 
x( p (x,Ql) = 0, 

xv p (x,Ql) = N v x a {l-xf, (47) 



xG p (x,Q 2 ) = N G xv p (x,Q 2 )=N G x a (l 
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where Nq = NqN v . Like in the LO analysis, [1-3], the input gluon distribution is 
proportional to the valence one and both have the valence-like form. Further, we 
impose two constraints given by Eqs. (J45)) and (|46|) in both types of models: FFNScxr- 
and CJK NLO. These constraints allow us to express the normalization factors N v and 
Nq as functions of a, (3 and k. This leaves these three parameters as the only free 
parameters to be fixed in the fits to the F^ experimental data. 



5.5 a s running and values of 



We distinguish between the number of active quarks in the photon, Nf, and the number 
of quarks contributing to the running of a s , which we denote by N q , [2] . The Nf depends 
only on the choice of the type of the model: ones we decide to use the FFNScjk or 



CJK NLO approach, then Nt = 3 or 5, respectively, through the whole Q 2 range of 



evolution. On contrary N q depends on the Q 2 value. It is equal 3 when Q 2 < m 2 
and increases by one unit whenever Q 2 reaches a heavy-quark threshold, i.e. when 
Q 2 = ml 

The running of the strong coupling constant in NLO is given then by the well-known 
formula: 

a { s Nq) (Q 2 ) = 1 _ Pi lnln(Q 2 /AW) 2 ) 

4tt " (3 \n(Q 2 /M N ^ 2 ) (31 (A^) 2 ) 2 1 ' 

with 

A) = ll-|iV ff , Pi = 102 - yiV 9 . (49) 

For N g = 4 we took the QCD scale equal to 280 MeV [22]. In order to 
ensure the continuity of the strong coupling constant at the heavy-quark thresholds the 
condition as\m\) = ai Nq+l \m 2 h ) is imposed. We calculate values of the remaining 
\( n q) constants according to the above condition. We utilizing the relation between 
constants A™ and A^ 1 ) given in [21]. That way we obtain A^ 3 ) = 0.323 MeV and 
A( 5 ) = 0.200 MeV. 

The parton evolution equations depend on N q through their dependence on a s (Q 2 ) 
and independently p and Pi, see Eqs. (|37jl and (|38|) . Therefore, because of the implicit 
introduction of the heavy-quark thresholds into the a s running we must proceed in three 
steps to perform the DGLAP evolution. In the first step, describing the evolution from 
the input scale Qo to the charm-quark mass m c , the hadronic input <2h ad (a;, Qo) ^ s taken 
from the VMD model. In the second step we evolve the parton distributions from m c 
to the beauty-quark mass, m^. Then, a new hadronic input is given by the sum of the 
already evolved hadronic and point-like contributions. The point-like distribution at 
Q 2 = m 2 becomes zero again. The same is repeated for Q 2 > m^. 



6 Global fits - results 

We have performed a series of fits to the F^x^Q 2 ) data [23]- [36]. All together 192 
data points were used, including the recent results of the ALEPH Collaboration [34], 
which replaced the old prelimiary data [38] , as well as recent preliminary results of the 
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DELPHI Collaboration [31] replacing the old prelimiary data [37]. In that case only 
the LEP1 data were applied because no final results for the LEP2 data are given in 
a form useful for our global fits 4 in [31]. Our fits based on the least square principle 
(minimum of \ 2 ) were done using Minuit [42]. Systematic and statistical errors of 
data points were added in quadrature. 

First, we made two test fits in which, following the analysis of the GRS group, [8], 
we assumed the input scale value Ql = 0.4 GeV 2 . The results of both fits are presented 
in table 01 The first two columns show the quality of the fits, i.e. the total y 2 f°r 192 
points and the y 2 P er degree of freedom. The fitted values for parameters a, (3 and k 
are presented in the middle of the table with the symmetric parabolic errors obtained 
from Migrad requiring Ay 2 = l 5 . In addition, the values for N v and N g obtained 
from these parameters, using the constraints (J45j) and (|46j) . are given in the last two 
columns. As can be seen the fits give high \ 2 per degree of freedom. Further, the 
relatively small k value indicates respectively small contribution of heavier mesons 
which is compensated by a very high, especially in case of the CJK NLO model, 
valence-quark input (N v = 9.26). 



NLO models 


x 2 


X 2 / DOF 


K 


a 


P 


N v 


N G 


FFNS CJi ,l 


318.0 


1.68 


1.40±0.06 


1.15±0.16 


1.24±0.42 


2.80 


0.949 


CJK 


299.9 


1.51 


1.40±0.05 


1.38±0.18 


3.40±0.77 


9.26 


2.19 



Table 3: The x 2 an d parameters of the test fits for 192 data points for FFNScjxI NLO and 
CJK NLO models with assumed Qq = 0.4 GeV 2 . The a, j3 and k symmetric parabolic errors 
obtained from Migrad requiring Ay 2 = 1. 



In order to obtain better agreement with the data two other fits were performed, in 
which Qq was treated at first as a free parameter. For both models we obtained very 
similar Q 2 « 0.765 GeV 2 (Q = 0.875 GeV) values. Therefore in our final CJK and 
FFNScxkNLO fits we fixed the input scale as Q 2 = 0.765 GeV 2 . The Q = 0.875 GeV 
value agrees with the Q = 0.85 ± 0.09 GeV obtained in the recent NLO fit presented 
in [10]. The results of both CJK and FFNS^xftT NLO fits are presented in table EJ 
The errors shown in the table were obtained from Minos requiring Ay 2 = 1 . 

As expected, estimation of the Qq value through the test fits, allowed us to find sets 
of parameters which give better agreement of both models with the data. We use these 
sets as the final results of our analysis. Comparison of the x 2 1 dof values obtained 
in the test and final fits shows the ~ 0.40 and ~ 0.30 improvment in the FFNScjkT 
and CJK NLO models respectively. Further, in both cases we observe an increase of 
the k parameter with the simultaneous decrease of and N v . That change can 
be explained by the fact that in higher Q 2 valence distribution broadens (change of 
«,/3) but at the same time the constraint (|45|) maintains its integral constant. Finally, 

4 All LEP2 data are given in three sets obtained with the TWOGAM [39], PHOJET [40] and 
PYTHIA [41] Monte Carlo generators. 

5 For test fits we applied Migrad instead of Minos which requires much more computer time 
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the Nq parameter grows at higher because the gluon contribution to the photon 
structure increases with the increasing Q 2 . 

The choice of the F^x^Q 2 ) distributions predicted by the test and final fits are 
presented in Figs. |U and El for the FFNScxrT and CJK NLO models, respectively. As 
the main difference we observe faster growth of the F^ function at small x in the case 
of the test fits with = 0.4 GeV 2 , especially at small Q 2 . The medium and high-x 
predictions for the photon structure function are similar apart from the low Q 2 region. 



NLO models 


x 2 


X 2 / DOF 


K 


a 


P 


N v 


N G 


FFNScjkI 


243.3 


1.29 


9 9CO+0.108 
z.zoo_ 096 


c i02 +0 - 071 


0.690±g;gg 


0.685 


2.369 


CJK 


256.8 


1.37 


2 662 +a108 


4Ofi+0.063 


1-0131^ 


0.745 


3.056 



Table 4: The x 2 an d parameters of the final fits for 192 data points for FFNScjftrl NLO 
and CJK NLO models with assumed Qq = 0.765 GeV 2 . The a, (3 and k errors obtained from 
Minos requiring A% 2 = 1. 



Moreover, we see that unlike in the LO case (see [1] and [2]), and in the case of the 
test fits with Ql = 0.4 GeV 2 , the obtained \ 2 P er degree of freedom is better in the 
standard type FFNS model than in the CJK model. The origin of the relatively high 
X 2 1 dof in both fits comparing to much lower \ 2 « 0.9 obtained in [10] is the same as 
in the LO case, discussed in detail in [3]. Namely, the exclusion of the CELLO [23] 
and recent DELPHI [31] data from the fits leads to much better agreement among fits 
and the data. 

Further, we performed test fits with the modified set of the input densities (I47|) . 
Namely, we used the input gluon distribution independent of the valence one 

x( p (x,Ql) = 0, 

xv p (x,Q 2 ) = N v x a {l-xf, (50) 
xG p (x,Q 2 ) = N G x aG {l-xY G . 

Fits gave very high values of the «g and j3 G parameters, around 14 and 16 respectively, 
leading to fast oscillations of the input gluon densities and all parton densities at low 
Q 2 . The Nq parameter calculated using constraint (|4l)|) is of the order 10 9 . Moreover, 
the «g and (3g values have very high uncertainties and therefore it is hard to obtain 
the final convergence of the fits. Taking all this and the fact that the x 2 1 dof values 
obtained in the test fits are slightly higher than the once presented above, we conclude 
that the choice of the simple form of the input distributions of Eqs. (}4T|) is very 
reasonable even at ~ 0.8 GeV 2 . 

6.1 Comparison of the CJK and FFNS CJir NLO fits with F 2 7 
data and other NLO parametrizations 

Figures EHU show a comparison of the CJK and FFNScjaT NLO fits to F^(x, Q 2 ) with 
the experimental data as a function of x, for different values of Q 2 . Also a comparison 
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with the GRS NLO [8] and AFG NLO [43] parametrizations is shown. (If a few values 
of Q 2 are displayed in a panel, the average of the smallest and biggest one was taken 
in the computation.) As can be seen in Figs. El and both CJK and FFNScxr-1 NLO 
models predict a much steeper behaviour of the F^x.Q 2 ) at small x with respect 
to other parametrizations. In the region of medium and high x, the behaviour of 
the F^x^Q 2 ) obtained from the FFNScxk-I and CJK NLO fits is similar to the one 
computed with the GRS NLO parametrization while it is very different to the AFG 
parametrization prediction, see Figs. |S1 and El Below the charm-quark threshold the 
CJK model gives the highest F^ and the GRS NLO parametrization the lowest F^ of 
the three similar predictions while above the threshold that order is inverted. Finally, 
at high Q 2 and high x the F 2 7 computed in the FFNScj^I NLO model and GRS NLO 
parametrization become undistinguishable while the CJK NLO prediction is visibly 
larger. 

Next, in Figs. EH and HU we compare the CJK and FFNS CJK 1 NLO fits with the 
corresponding CJK and FFNSc*jk2 LO fits presented in [2] . Figure El shows very 
small differences, only in the low x region the FFNScj/fl NLO model gives greater 
values then other models and parametrizations. On the other hand, as it can be seen 
in Fig. ^2 the NLO models predict higher F% than the LO ones at the x region right 
below the charm-quark threshold. Moreover the NLO predictions vanish more rapidly 
at x — > 1. 

Apart from the direct comparison with the photon structure-function data, we 
perform another comparison, this time with LEP data that were not used directly 
in our analysis. Figures IT21 and [T31 present the predictions for i^x, Q 2 ), averaged 
over various x regions, compared with the recent OPAL data [36]. For comparison, the 
results from the GRS and AFG NLO parametrizations are shown as well. In Fig. IT2*lwc 
observe that in the medium-x range, 0.1 < x < 0.6, all models give expectations of the 
same shape, all in fairly good agreement with the experimental data. Still, in small Q 2 
the Fl computed with the AFG NLO parametrization is above the others, and starting 
at approximately Q 2 = 3 GeV 2 the predictions of the GRS NLO parametrization are 
lower than the other three. In the x ranges presented in Fig. ^]we also notice many 
differences among the averaged F^(x, Q 2 ) predictions computed in various models. The 
only exception give the FFNScjjfl NLO and the GRS NLO predictions which are very 
close to each other in the whole Q 2 range and in all, apart from the lowest, x ranges 
examined. Let us now name all the dissimilarities observed in FiglTBl Firstly, the 
F^Q 2 ) predicted by the CJK NLO model shows clear change of the behaviour in the 
high Q 2 and high x range where we observe its larger increase than predicted by the 
other models. Secondly, the expectations of the AFG parametrization in the same 
regions (also for small Q 2 ) lie much below the other results. Finally, the GRS NLO 
parametrization predicts slightly lower F^iQ 2 ) values than the remaining models. 

6.2 Parton densities 

In this section we present the parton densities obtained from the CJK and FFNScj^I 
NLO fits and compare them with the corresponding distributions of the GRV [7] , GRS 
and AFG NLO parametrizations. First, we present all parton densities at Q 2 = 10 



21 



GeV 2 (Fig. ITljl . The biggest difference between our CJK NLO model and others is 
observed, as expected, for the heavy-quark distributions. Unlike for the GRV and AFG 
NLO parametrizations (there is no beauty-quark density in the AFG NLO parametriza- 
tion), our c 7 (x, Q 2 ) and 6 7 (x, Q 2 ) densities vanish not at x — 1 but, as it should be, at 
the kinematical threshold. Moreover, we notice a huge difference between the beauty- 
quark distributions computed in two models. The GRV NLO parametrization predicts 
the 6 7 (x, Q 2 ) values of the same order as the c 7 (x, Q 2 ). In Fig. El where the charm- 
quark distributions are presented for various Q 2 values, apart from different threshold 
behaviour we observe that CJK c 7 (x, Q 2 ) distribution is the largest of all down to 
very small x value (depending on Q 2 ) where the GRV density becomes larger. The 
light-quark and gluon densities calculated in our models and the GRV and GRS NLO 
parametrizations generally have the same shapes but their values differ. Main differ- 
ences are observed in the up-, down-quark and gluon distributions. The u 1 densities 
computed in both our models are lower than in other parametrizations at high x. In 
the Fig. Elwe see differences among predictions of all models at very low x. Still we do 
not find a general pattern of those distinctions. In the case of the gluon distribution, 
see Figs. El and El we observe that at all Q 2 regions the CJK model gives the largest 
predictions. At high Q 2 , G 7 (x, Q 2 ) calculated in the FFNScxr-1 NLO model and in 
GRV NLO parametrization become similar. The GRS and AFG NLO predictions lie 
below all other lines in the whole x and Q 2 range. 

Further, we compare our CJK and FFNScjifl NLO parton distributions with the 
corresponding CJK and FFNScxr-2 LO ones. As can be seen in Fig. El for Q 2 = 10 
GeV 2 , the NLO models predict higher up- and down-quark densities in the medium- a; 
region and steeper increase of the gluon distribution. The growth of G 7 (x, Q 2 ) at low 
x is especially fast in the CJK NLO model case while the FFNScjftrl NLO model 
prediction tends to fuse with the LO ones. Moreover, all NLO quark densities vanish 
more rapidly at x — > 1 than the LO quark densities. 



6.3 Comparison with R 



7 

2,c 



In Fig. El we present the CJK NLO predictions for the F 2 7 c . For Q 2 = 5,20, 100 and 
1000 GeV 2 we compare the individual contributions included in the model. Analogically 
to the LO case, see [1] and [2], we introduce the following notation 



2xe 2 c c-<(x,Q 2 )+FZ c \ dil (x,Q 2 



2,clres 



X, Q ) -f2,clres,sub(-£> Q )" 



2c * 



(51) 
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with 



^ c k(*,Q 2 ) = 3e 2 " 



7T 



Am 2 ( M o ^Am 2 c 2 % m i\ , 1 + 



F-^U^Q 2 ) = f -G-'(y,Q 2 )CiS } (-,%), (52) 

Jxc y y m c 

y 



^ c | rc s,sub(x, g 2 ) = In g • 2 /* ^(y, g 2 )P f g 



(*,g 2 ) = / ^(c^ + c^^g 2 )^^). 

y y 



Let us stress that the above direct contribution effectively includes both the direct 
and direct subtraction terms, as discussed in detail in section l4~4l Therefore, it is not 
the same as the .F^cldir^j Q 2 ) in the LO analyses. That way in the NLO analysis we 

removed the troublesome subtraction term ~ . Unfortunately, also the 2c * C^ q 
contribution does not vanish at the threshold and therefore an additional constraint 
2c * C^l = for Xc > 1 is necessary. This term is important in the high and medium 
x, its values in the latter range are negative. As can be seen in the plot all other 
terms vanish in the W — *> 2m c threshold. The direct term influences the F 2c only in 
the midium and high x at low-g 2 range, where it is negaive. The charm-quark density 
contribution, i.e. the term 2xe 2 c y (x, Q 2 ), dominates the F 2c in the whole kinematically 
available x range. In the low-x region also the resolved-photon contributions increase, 
but they cancel each other. 

A good test of the charm-quark contributions is provided by the OPAL measurement 
of the F^ C i obtained from the inclusive production of D*^ mesons in the photon-photon 
collisions [44]. The averaged F^ c has been determined in the two x bins. These data 
points are compared to the predictions of the CJK and FFNScjkI NLO models, as 
well as of the corresponding CJK and FFNScje'2 LO models in Fig. EH The CJK 
NLO model seems to give the best description of the data, especially for the high-x 
bin. The FFNS type models predict the same F 2c at x > 0.03, while below that value 
the FFNScxr-2 LO predicts higher F% c which finally becomes very similar to the NLO 
results. 



6.4 Test of importance of 0(a 2 s ) and 0(aa s ) terms 

In this section we would like to present the results of the fit of the FFNScja-2 NLO 
model to the data. As discussed in Sections 13.21 and 13.41 in this approach we include 
the 0(af) and 0(aa s ) contributions to the photon structure function: 

1 

x 



-F 2 >,g 2 




(53) 



E(^ + 97)(2/>Q 2 )(^c 2 ( ; 



V m h 



2m \ i I)) 

y ™ h 



+e 2 h G^y,Q 2 )C h 2 S\ X -^] 

y " L h 



, aa s (Q 2 ) 4^,(1), 



'2,7 



X, 
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Alike in the main CJK and FFNScjx 1 NLO models our fit was performed in two 
steps. First, we treated Ql as a free parameter in order to obtain better agreement 
with the data, next it was fixed to the estimated value, Ql = 0.716 GeV 2 , and a final 
fit was done. Its results are presented in table all errors were obtained from Minos 
requiring A% 2 = 1. 



NLO model 


x 2 


X 2 1 ' DOF 


K 


a 


P 


N v 


N G 


FFNS C j^2 


252.2 


1.31 


1 892 +0 ' 101 
i.oyz_ 090 


n c-97+0.089 
U.OZrf _ .080 


0.966lgS 


0.797 


2.727 



Table 5: The x 2 an d parameters of the fit for the FFNScjk2 NLO model with assumed 
Qq = 0.716 GeV 2 . The a, (3 and tz symmetric parabolic errors obtained from Minos requiring 
A X 2 = 1. 



First let us compare the results of the two FFNScxr: NLO models, given in ta- 
bles |U and El respectively. We notice that higher x 2 1 dof is obtained in the model 
including additional terms. Next, the input scale and the parameters obtained in both 
approaches are very similar. Largest distinction is observed between the k values and 
can be explained by the small difference between the Ql values (0.716 and 0.786 GeV 2 , 
respectively). Namely, when we start the evolution of the parton densities at lower 
input scale then smaller input distributions are required in order to obtain the same 
results at the given Q 2 . 

Further, we compare results computed using four various FFNScja: models, the 
NLO ones described in this work as well as the FFNScjkI & 2 LO models, presented 
in [2]. The main differences among the four models, apart from the different order of 
the QCD DGLAP equations applied, lie in the heavy-quark contributions to F^. Those 
differences are summerized in table|H]where all included or discluded terms proportional 
to Wilson-coefficient functions are listed. 



FFNS C jk model 


n h,(o) 

°2,7 






L/ 2,G 




U 2,q 


LO 1 


+ 












LO 2 


+ 


+ 










NLO 1 


+ 


+ 










NLO 2 


+ 


+ 


+ 


+ 


+ 


+ 



Table 6: The heavy-quark terms proportional to the Wilson coefficient functions contributing 
(or not) to the heavy-quark structure function in various FFNScj/f models. 

In Fig. we compare the parton densities computed in the FFNScjx models 
obtained in LO and NLO. We observe that two LO models predict nearly the same 
parton distributions which begin to differ in the very low-x region (the larger differences 
at x below 0.1 are not shown). The distinctions between the pairs of LO and NLO lines 
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are large at high x where, as discussed before, the NLO distributions decrease much 
earlier. The maximal values of the NLO u 1 and s 7 distributions are reached at smaller 
x, x ~ 0.75 than in the LO case where the corresponding x value is x ~ 0.95. Moreover, 
there exist important differences between the two NLO fits, observed especially in the 
up- and down-quark distributions, which indicate the importance of the higher order 
0(a 2 ) and 0(aa s ) contributions to F^(x, Q 2 )&t medium x. Smaller difference between 
the LO results seen in Fig. |!H]may be explained by the fact that in that case the models 
differ only by the resolved-type ~ C\ 'jp contribution, which plays important role only 
in the low-x region. Very similar conclusion may be drawn from the comparison of the 
F^{x, Q 2 ) distributions plotted in Figs. 12*2*1 and 12*31 Alike in the parton densities case 
the LO F^x^Q 2 ) distributions differ only in the very low-x region, see Fig. "211 while 
the NLO ones mostly at high x. 

6.5 Influence of the DELPHI LEP2 data 

Finally we tested the possible influence of the recent DELPHI data [31], which were 
not included in the final fits, on the results of our models. We performed test fits for 
the CJK NLO model including separately each of the sets of the LEP2 data obtained 
with TWOGAM [39], PHOJET [40] and PYTHIA [41] Monte Carlo generators and 
presented in [31]. As can be seen in Figs. "2*4*1 and "2*51 for the choice of the F^x^Q 2 ) 
distributions and the parton densities computed at Q 2 = 10 GeV 2 , respectively, the 
results of that fits are very similar to the expectations of the final CJK NLO fit (denoted 
as "no LEP2"). Still, the \ 2 1 dof obtained in the test fits are much worse, see tabled 



set 


x 2 


# of points 


X 2 / DOF 


TWOGAM 


307.4 


209 


1.50 


PHOJET 


312.0 


207 


1.54 


PYTHIA 


341.1 


209 


1.66 



Table 7: Number of experimental points used in the test fits for the CJK NLO model and 
the resulting \ 2 - 



7 Disscusion and conclusions 

A new analysis of the radiatively generated parton distributions in the real photon 
based on the NLO DGLAP equations is presented. An updated set of the F^x^Q 2 ) 
data including recent ALEPH as well as new preliminary DELPHI measurements has 
been used to perform three global fits. Our models base on two schemes, the Fixed 
Flavour-Number Scheme and Variable Flavour- Number Scheme (ACOT(%)) which 
were already applied in our former LO analysis, see [1]- [3]. We observe that in com- 
parison to the LO models higher input scale of the NLO DGLAP evolution is required 
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in order to satisfactory describe the data. Further, unlike in the LO case we obtained 
lower x 2 / dof m the FFNS type approach than in the CJK model. The test fits showed 
that a change of the simple form of the input parton distributions applied in [1-3] and 
in this work does not lead to the improvement of the fits. 

Predictions of the FFNScxr-1 and CJK NLO fits are compared with the F%(x, Q 2 ) 
data, other parametrizations and corresponding LO models. Both models describe 
very well the Q 2 evolution of the F^(x, Q 2 ) , averaged over various x-regions. The F^ c 
calculated in the CJK NLO model gives best agreement with the data. 

Next, we check the difference between the results obtained in the FFNScj/fl NLO 
model and the FFNScxfi:2 NLO model which includes additional 0(a 2 ) and 0(aa s ) 
contributions to F^(x, Q 2 ). Those higher order terms prove to be of importance. There- 
fore, further studies of the CJK model including terms of the same order are required. 

Finally, we examine the influence on the results of our models of the various sets of 
the LEP2 data presented in the recent preliminary results of the DELPHI Collaboration 
but discluded from final fits. Test fits show very small change of the F^x^Q 2 ) and 
parton distributions obtained using these data comparing to the final fits but large 
deterioration of their quality. 

Fortran parametrization programs for our CJK and FFNS^jif NLO models can be 
obtained from the web-page [45]. 
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Appendix: The Wilson-coefficient functions 

Lists of all Wilson-coefficient functions used in this analysis are given in tables ^ and El 
The Wilson coefficients C^q, C^G and C^\L have been first calculated in the operator 
product expansion [46] in the Mellin space in the articles [47] and [48] , respectively. In 
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the (x, Q 2 ) space the C2I in the MS factorization scheme reads (see for instance [8]) 



Cg(*)=6 



-1 + 8x(l -x) + (x 2 + (1 - x) 2 ) In 



1 — x 



X 



(54) 



The lowest order heavy-quark coefficient C^'^, given by the Bethe-Heitler 7* + 7 



q + q process, has been first presented in [49] . It has the following form: 



r ,h,(o)f Q 2 > 

( - y 2,7 l X ' 2' 
m h 



7T 



1 + 8x(l — x) — x(l — x) 



(55) 



+ ln 



l + (3 
1-P 



Am? 

x 2 + (l-x) 2 + x(l-3x)—± 

v 



x 



imj 



Q 4 



-, 4m? x 
1 — 



with /3-yi- (jz^. 

When we neglect the heavy-quark mass in some parts of the coefficient C^ '' we 
can rewrite it as in Eq. (J21j) 



a 



?i,(0),MS/„ Q 2 



Cg' MS (a;) + e^)' MS (a;)ln 



Q 2 



mt 



(56) 



+ 6e 4 h 



4m 



Q 



Ami 



-(3x(l - x)—± + x(l - 3a;)— f - x 



Q 7 



Q 4 



in 



1-/5 



where the direct photon-quark splitting function is 



= 3[x 2 + (1 - x) 



(57) 



The q 1 function related to the 7* + G — > g + g process reads 



/?)7 



6 



2 " V 'm 2 , 



(58) 



a 



7T 

+ ln 



1 + Sx(l — x) — x(l — x) 



Am 2 



Q 2 



Am 

x 2 + (l-x) 2 + x(l-3x)—± 

v 



x 



iml 



Finally, the higher order heavy-quark coefficients, C^'i , C^Gj ^2,5 an d are 
calculated in the serious of publications [13]. They are given in the form of a two 
dimensional array. Those 0(a 2 ) and 0(aa s ) order Wilson coefficients are given as 



y(2) 



2,i v X ' 2 / 2 ' 



(59) 



where /1 2 is the mass factorization scale. In our calculations we simplify the above 
formula by choosing fi 2 = m\. 
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Figure 4: Predictions for the F^x, Q 2 )/a for the FFNScxkI NLO model obtained by 
the main fit with Q 2 = 0.765 GeV 2 and by the test fit with Q 2 = 0.4 GeV 2 . 
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Figure 6: Predictions for the F 2 7 (x, Q 2 ) /a for the CJK NLO and FFNS C j^ NLO models 
and GRS NLO [8] and AFG NLO [43] compared with the experimental data [23]- [36], 
for small and medium Q 2 as a function of x (logarithmic scale). If a few values of Q 2 
are displayed in the panel, the average of the smallest and biggest Q 2 was taken in the 
computation. or, 
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Figure 7: The same as in Fig. El for Q 2 > 20GeV 2 
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Figure 9: The same as in Fig. [7|for a linear scale in x. 
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Figure 10: Predictions for the F^{x,Q 2 )/a for the CJK NLO and FFNS C j^ NLO 
models compared with corresponding LO fits from [2] and with the experimental data 
[23]- [36], for small and medium Q 2 as a function of x (logarithmic scale). If a few 
values of Q 2 are displayed in the panel, the average of the smallest and biggest Q 2 was 
taken in the computation. w 
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Figure 11: The same as in Fig. E3for a linear scale in x. 
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Figure 12: Comparison of the OPAL data [36] for the <3 2 -dependence of the averaged 
over 0.1 < x < 0.6 F^/a with the predictions of the CJK NLO and FFNS C jk NLO 
models. 
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Figure 13: As in Fig. IT^l for F^x^Q 2 ) /a, averaged over four different x ranges. 
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Figure 14: Comparison of the NLO parton densities predicted by CJK NLO and 
FFNScxkI NLO models and by GRV NLO [7], GRS NLO [8] and AFG NLO [43] 
parametrizations at Q 2 = 10 GeV 2 , function of x. 



39 



10 3 



CJK NLO 




Figure 15: Comparison of the up-quark density at four values of Q 2 in the CJK and 
FFNS C j^l NLO models with the GRV NLO [7], GRS NLO [8] and AFG NLO [43] 
densities. 
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Figure 17: The same as in Fig. 1151 for the charm-quark density. 
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Figure 18: Comparison of the parton densities predicted by the CJK NLO and 
FFNScxr-1 NLO models and by corresponding LO models, CJK LO and FFNScxr-2 
LO (i.e. including the resolved-photon j*G — > hh contribution) from [2] at Q 2 = 10 
GeV 2 , as a function of x. 



43 




Figure 19: Comparison of various contributions to the photon structure function 
F^ c (x, Q 2 )/a in the CJK NLO model for Q 2 = 5, 20, 100 and 1000 GeV 2 . 
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Figure 20: Comparison of the structure function F2 C (x,Q 2 )/a calculated in the CJK 
and FFNSc-jk models for NLO & LO with the OPAL measurement [44]. 
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Figure 21: Comparison of the parton densities predicted by the FFNScxk-1 & 2 NLO 
models and by the corresponding FFNScj^I & 2 LO models presented in [2] at Q 2 = 10 
GeV 2 , as a function of x. 



15 



0.4 
0.2 




Tnj 1 rTTTTTTj 1 rTTTTTTj r 


1 1 1 llll 


: i lo - 




r FFNS C j^ 2 LO - 




: 1NLO-- 




: 2 NLO - 




mi i i i 






0.8 p 
0.6 r 
0.4 ^ 
0.2 t 


0.8 



0.6 - 
0.4 r 
0.2 t 




I I llllj 1 TTTTTTTj - 

OPAL 3.7 



i 1 1 iiiij 1 i i i niij r 

AMY 6.8 
CELLO 7.0 ■♦■ 



1 1 niij — i — i 1 1 niij 

TPC 1.3 




1 1 1 nij — i 1 1 1 1 nij — 

L3 1.9 
OPAL 1.9 ■♦■ 




1 1 1 niij — i — 1 1 1 niij 

PLUTO 4.3 



TTlj 1 I I I I lllj 1 TTTTTTTj 1 - 

L3 5.0 
TPC 5.1 ♦■ 
TOPAZ 5.1 ^ 
DELPHI 5.2 O 




ti] — i — i 1 1 niij — i — i 1 1 niij — 

OPAL 8.9 
OPAL 9.0 ■♦■ 
PLUTO 9.2 t 




TTTj 1 I I I I lllj 1 I I I I lllj 1 - 

ALEPH 9.9 
i- OPAL 10.7 ■♦■ 

L3 10.8 




1 


ji.j , , 


1 1 niij — i i i i niij — r 




TTTI 


0.8 




DELPHI 12.0 


▲ 






DELPHI 12.7 


■♦■ 




0.6 


T 

:\ 

i— V 


CELLO 13.1 


T 




0.4 




^V* TAJ 
















0.2 




▲ 









~;..\ . . 


1 1 i_ 







TTTTTTTj 1 — I I I llllj 1 - 

ALEPH 17.3 -± 
OPAL 17.5 ♦ 
OPAL 17.8 t 




jm 1 i 1 1 1 1 mi 1 i i 1 1 1 hi 1 — nn 

\ TASSO 23.0 ^ ; 
-\ , L3 23.1 id 
JADE 24.0 TlV 




2 
1.5 
1 

0.5 




ij 1 i 1 1 niij 1 i i i niij n 

OPAL 59.0 -a- 
ALEPH 67.2 -♦- 




0.001 0.01 0.1 




■nn — i — i 1 1 niij 1 — i 1 1 niij r- 

:'v\ AMY 390 
-_\\ DELPHI 400 ♦ 




uL 



uL 



uL 



0.001 0.01 



0.001 0.01 0.1 



X 



Figure 22: Predictions for the F^{x, Q 2 )/oi for the FFNScjit models, as in Fig. EU 
compared with the corresponding LO fits from [2] and with the experimental data [23]- 
[36], for small and medium Q 2 as a function of x (logarithmic scale). If a few values of 
Q 2 are displayed in the panel, the average of the smallest and biggest Q 2 was taken in 
the computation. 
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Figure 23: The same as in Fig. |22]for a linear scale in x. 
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Figure 24: Comparison of the parton densities predicted by the test fits obtained with 
the CJK NLO model using various sets of the DELPHI LEP2 data [31] at Q 2 = 10 
GeV 2 , as a function of x. 
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Figure 25: Predictions for the F^x^Q 2 ) / a for the CJK NLO model obtained in the 
test fits using various sets of the DELPHI LEP2 data [31]. 
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